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A systematic method to calculate the low energy 
spectrum of 5(7(2) Yang-Mills quantum mechanics 
with high precision is given and applied to obtain the 
energies of the groundstate and the first few excited 
states. 
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The so-called Yang-Mills mechanics originates 
from Yang-Mills field theory under the supposi- 
tion of the spatial homogeneity of the gauge fields. 
Since its introduction twenty years ago [jj] this 
model has been studied extensively from different 
points of view § - The model is known to ex- 
hibit stochastical behaviour at the classical level. 
In particular it has been shown for a simplified two- 
dimensional version of the model that it is nonin- 
tegrable [||. How this reflects itself on the corre- 
sponding quantum mechanical problem is still un- 
der investigation |Q - [gj , , @ • Furthermore it 
offers the possibility to develop methods gcneraliz- 
able to the study of the corresponding non-Abelian 
gauge field theory. With this point of view Yang- 
Mills quantum mechanics has been investigated in 
the small volume approximation in the small 
coupling limit on the lattice using the semiclassi- 
cal approximation Q , and to estimate the glueball 
spectrum [ fhlf . Furthermore the question of the ex- 
istence of normalizable zero-energy wave functions 
in supersymmetric matrix models Jll| as general- 
izations of Yang-Mills mechanics is under intensive 
investigation (see e.g. |l7| and references therein). 

In the present work we would like to investigate 
the groundstate and the first few excited states of 
SU(2) Yang-Mills quantum mechanics using the 
variational approach. The quantum treatment of 
such models is confronted with the difficulties due 
to the existence of the non-Abelian Gauss law con- 
straints and the presence of the unphysical degrees 



of freedom. In jl3| it has been shown that classical 
SU(2) Yang-Mills mechanics, after the elimination 
of all unphysical degrees of freedom, reduces to the 
equivalent unconstrained system described by the 
Hamiltonian 
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This Hamiltonian is defined on the phase space 
spanned by the three canonical pairs Xi,Pi and 
the three 5*0(3) left-invariant Killing vectors & 
satisfying the Poisson bracket relations , £j- } = 
—^ijk^k ■ Each of the three configuration space 
variables Xi is defined on the positive half line. The 
Hamiltonian (|]) has been obtained by writing the 
spatially homogeneous SU(2) gauge fields Af(t) in 
the so-called polar representation Jl8| 

A-ai (q, Q) = Oak (q) Qki , 

with the orthogonal 3x3 matrix O parametrized 
by the three angles q and the positive definite 3x3 
symmetric matrix Q. The three angles q are pure 
gauge degress of freedom, which do not enter the 
Hamiltonian while their conjugated momenta are 
vanishing according to the Gauss law constraints. 
The variables Xi appearing in the Hamiltonian ([!]) 
are the eigenvalues of the positive definite symmet- 
ric matrix Q 

Q = R T ( X )D(x 1 ,x 2 ,x 3 )R(x) ■ 

The angles x m the orthogonal matrix R(x) and 
their canonical momenta are combined in the 
Killing vectors . We emphasize that by defi- 
nition of the polar decomposition the diagonal el- 
ements Xi are positive definite which has impor- 
tant consequences on the quantum level, as will be 
shown below. 

The transition to the corresponding quantum 
Hamiltonian is free of operator ordering ambigu- 
ities. The operators Ji = Rij^j corresponding to 
the right-invariant 50(3) Killing vectors commute 
with the Hamilton operator. Its eigenstates can 
therefore be classified according to the quantum 
numbers J and M as eigenvalues of the total spin 
f 2 and its projection J3. We shall consider here 
only the spin zero sector, J = £ = 0, for which 
the Schrodinger eigenvalue problem reduces to 
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For the Hermiticity of the Hamiltonian Hq dehned 
in the wave functions have to satisfy the bound- 
ary conditions 



lim , t>(x 1 ,x 2 ,x 3 ) = , 
lim \Jj(xi,X2,x 3 ) = 0. 
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The potential in (Q) has three flat valleys of zero 
energy, x\ = x 2 = and x% arbitrary, and the two 
others obtained from this by cyclic permutation. 
Close to the bottom of the valleys the potential is 
that of a harmonic oscillator with a width narrow- 
ing down for larger values of X3. 

One of the main results known for such a prob- 
lem is that it has been proven to have a discrete 
spectrum due to quantum fluctuations g, al- 
though the classical problem allows for scattering 
trajectories. Based on the well-known operator in- 
equality [] 
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it follows that 



H > - (-A + 3V2 5 (xi + x 2 



•2:3. 



=: -H' 
2 



(5) 



Since the Hamiltonian H' is known j|] to have 
a discret spectrum, this is true also for Ho- An 
important open question is at which energy the 
groundstate is. The well-known nonnormalizable 
zero energy solution ^e=o '■— A ex p[ — 9x1X2X3} 
p9[ , p0| , does not satisfy the necessary bound- 
ary condition (||). The knowledge of the ground- 
state energy of H' in inequality (J5J) would provide 
a lower bound for the groundstate energy of Hq. 
Due to the additive structure of the potential term 
in H' one can make a separable ansatz for the solu- 
tion of the corresponding eigenvalue problem. The 
energy of the lowest such separable H' eigenstate 
satisfying the above boundary conditions (0) and 



E'=6\(o\^g/2) 2/3 = 9.1924 
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where £ — —2.3381 is the first zero of the Airy 
function. We shall see in this letter that there exist 



1 The number 3 in r.h.s of the inequality originates 
from the positiveness of x. 



nonseparable solutions which have an energy below 

To obtain an upper bound £ for the groundstate 
energy Eo of the Hamiltonian the most powerful 
tool is the Rayleigh-Ritz variational technique plj , 
based on the minimization of energy functional 
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Guided by the harmonic oscillator form of the val- 
leys of the potential in (]|) close to there bottom 
a simple first choice for a trial function compat- 
ible with the boundary conditions (^) and (||) is 
to use the lowest state of three harmonic quantum 
oscilators on the positive half line 
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The stationarity conditions for the energy func- 
tional of this state, 

OFT 1 V** (?> 9 2 1 

^ V 4 8 WjWk 

cyclic J 

lead to the isotropic optimal choice 

uj := cji = uj 2 = u> 3 = 3 V V /3 . (9) 

As a first upper bound for the groundstate energy 
of the Hamiltonian we therefore find 
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3 l/3 5 2/3 = 4 8676 ff 2/3^ 



This upper bound is in agreement with the lower 
bound of the energy functional for separable func- 
tions 
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obtained from the operator inequality (|5|) and the 
lower bound (||) for separable solutions of H'. 

In order to improve the upper bound for the 
groundstate energy of the Hamiltonian H we ex- 
tend the space of trial functions (@) and consider 
the Fock space of the orthonormal set of the prod- 
ucts 



(11) 



of the odd eigenfunctions of the harmonic oscillator 

( W /tt)V4 
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with the frequency fixed by (||). Furthermore the 
variational procedure becomes much more effec- 
tive, if the space of trial functions is decomposed 
into the irreducible representations of the residual 
discrete symmetries of the Hamiltonian (||). It is 
invariant under arbitrary permutations of any two 
of the variables PijXi = XjPij, PijPi = PjPij and 
under time reflections Txi = x{F, Tpi — —p{F , 



[H , Pi, 



0, 



[H ,T] =0. 



We shall represent these by the permutation oper- 
ator P12, the cyclic permutation operator P123 and 
the time reflection operator T, whose action on the 
states is 

P^^^i, 2^2,2:3) = ^(x 2 ,x 3 ,Xi) , 
Pi 2 ^(xi,x 2 ,x 3 ) = ^(x 2 ,xi,x 3 ) , 
T^{x 1 ,x 2 ,x 3 ) = ^*(xi,x 2 ,x 3 ) , 

and decompose the Fock space spanned by the 
functions ( |i"l| ) into the irreducible representations 
of the permutation group S3 and time reflection T. 
For given (rti, n 2 , n 3 ) we define 

— m 

* nnn ' — * nun > 

if all three indices are equal (type I), the three 
states (to = —1, 0, 1) 

1 2 
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when two indices are equal (type II), and the two 
sets of three states (m — —1,0,1) 

\U(™)± .— 
* n 1 n 2 n 3 • 

1 £ e -*W {Pl23 f (1 ± Pl2 ) *„ in2 „ 3 , 

* k=0 

if all (m, n 2 , n 3 ) are different (type III). In this new 

orthonormal set of irreducible basis states 4 , ^ rl ' )Q , 
the Fock representation of the Hamiltonian Hq 
reads 



The basis 
and P 12 T 
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The basis states \&jij^ a are eigenfunctions of Pi 
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Pl2T* = ±*S' V:C • ( 12 ) 

Under P\ 2 and T separately, however, they trans- 
form into each other 
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We therefore have the following irreducible repre- 
sentations. The singlet states \I/( ) + , the "axial" 
singlet states * (0) ~, the doublets 
and the "axial" doublets (*(+!)-; ^C" 1 )-). Since 
the partner states of the doublets transform into 
each other under the symmetry operations P\ 2 or 
T, the corresponding values of the energy func- 
tional are equal. 

Due to this decomposition of the Fock space 
into the irreducible sectors, the variational ap- 
proach allows us to give upper bounds for states 
in each sector. The values of the energy functional 
for the states in each irreducible sector with the 
smallest number of knots £[*j,oo + ] = 4.8676 g 2 / 3 , 
£[^ )+ ] = 7.1915 g 2 / 3 , 8[¥ %-\ = 13.8817 ff 2 / 3 , 
and £[tf£i2 1)_ ] = 15.6845 g 2 ' 3 give first upper 
bounds for the lowest energy eigenvalues of the sin- 
glet, the doublet, the axial singlet, and the axial 
doublet states. 

In order to improve the upper bounds for each 
irreducible sector, we truncate the Fock space at 
a certain number of knots of the wave functions 
and search for the corresponding states in the 
truncated space with the lowest value of the en- 
ergy functional. We achieve this by diagonalizing 
the corresponding truncated Hamiltonian fftmnk 
to find its eigenvalues and eigenstates. Due to the 
orthogonality of the truncated space to the remain- 
ing part of Fock space the value of the energy func- 
tional (Q) for the eigenvectors of -fftrunk coincides 
with the fftrunk eigenvalues. 

Including all states in the singlet sector with up 
to 5 knots we find rapid convergence to the follow- 
ing energy expectation values for the three lowest 
states Si, S 2 , S 3 

£[Si] = 4.8067 5 2/3 (4.8070 g 2/3 ), 
£[S 2 ] = 8.2515 5 2/3 (8.2639 g 2/3 ), 
£[S 3 ] = 9.5735 5 2/3 (9.6298 g 2/3 ), (13) 

where the number in brackets show the corre- 
sponding result for 4 knots. The lowest state Si is 
very close to the state ^qoo +: 



Si = 0.994562 * 
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0.0970056 



0.0145832 



(14) 



Similarly including all states in the doublet sec- 
tor with up to 6(5) knots the following energy 
expectation values for the three lowest states 
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£[Z^ ±1} ] = 7.1682 g 2/3 (7.1689 g 2/3 ), 
£[D ( 2 ±1] ] = 9.6171 g 2/3 (9.6394 g 2 ' 3 ), 
£[D 3 ±1] } = 10.9903 g 2/3 (10.9951 g 2/3 ). (15) 

have been obtained. Including all states in the ax- 
ial singlet sector with up to 8(7) knots we find the 
following energy expectation values for the three 
lowest states A\,A 2 , A3 

£[A X ] = 13.2235 g 2/3 (13.2275 g 2/3 ), 
£[A 2 ] = 16.6652 g 2 ' 3 (16.7333 g 2/3 ), 
£[A 3 ] = 19.1470 g 2/3 (19.3028 g 2/3 ). (16) 

Finally taking into account all states in the axial 
doublet sector with up to 8(7) knots we find the 
following energy expectation values for the three 
lowest states C^ l \c^ l) ,C^ l) 



,2/3x 



£ [C[ ±1} ] = 14.8768 g 2 ' 3 (14.8796 
£[C { 2 ±1} ] = 17.6648 g 2/3 (17.6839 g 2/3 ), 



£[C^] 



19.9019 g 2/3 (19.9914 
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We therefore obtain rather good estimates for the 
energies of the lowest states in the spin-0 sector. 
Extending to higher and higher numbers of knots 
in each sector we should be able to obtain the low 
energy spectrum in the spin-zero sector to high 
numerical accuracy. 

In summary comparing our results for the first 
few states in all sectors, we find that the low- 
est state appears in the singlet sector with en- 
ergy 4.8067 g 2 / 3 with expected accuracy up to 
three digits after the dot. For comparison we re- 
mark that due to our boundary condition (||) all 
our spin-0 states correspond to the 0~ sector in 
the work of jl0| where a different gauge invariant 
representation of Yang-Mills mechanics has been 
used. Their state of lowest energy in this sector 
is 9.52 g 2 / 3 . Furthermore in [H, using an anal- 
ogy of SU (N) Yang-Mills quntum mechanics in the 
large N limit to membrane theory, obtain the en- 
ergy values 6.4690 g 2/3 and 19.8253 g 2/3 for the 
groundstate and the first excited state. 

Finally we remark that an analogous calcula- 
tion shows that also the groundstate energy of the 
Hamiltonian H' in Eq. (||) is lower than the value 
E' sep of (||) for the lowest separable solution. 
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